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Abstract 



O ' In symmetrical nuclear matter the solutions of pion dispersion equation are investigated 

^ ! in the complex plane of the pion frequency uj. There are three well-known branches of 

\ solutions on the physical sheet : sound, pion and isobar, at the matter density less than 

i-j \ the critical one p < pc- At p > pc the fourth branch appears on the physical sheet. For 

this branch the condition cu^ < (in general case Rea;^ < 0) takes place. This points out 
the instability of the ground state which is possibly related to the pion condensation. 



in 



^ ■ PACS: 13.75.Gx Pion-baryon interactions; 21.60.Jz Hartree Fock and random-phase 

•jH . approximations; 21.65.+f Nuclear matter 

^ ! Keywords: pion dispersion equation, pion condensation, symmetrical nuclear matter 
O 

O ■ 1 Introduction 



Some time ago a bright phenomenon of pion condensation predicted in attracted the 
common attention and was widely investigated The reason for this prediction 

Ch ■ was the fact that one of the solutions of the pion dispersion equation in medium (let 

it be called Uc) turns to zero uj'^ik) = (at some momentum k ^ 0) while the density 
of medium increases 111]. This meant an appearance of excitations with zero energy in 



X 



' medium. In this case the ground state should be reconstructed in correspondence with 

^ phase transition. 

One variant of reconstruction was to include the pion condensate (taken in one or 
another form) into the ground state. This resulted in a prediction for pion condensation. A 
number of efforts was devoted to the search for this phenomenon. The result of discussions 
given in the book of A.B.Migdal et. al. p[ was that the pion condensation manifested itself 
weakly and probably was absent (at least not observed) at normal density. Nevertheless, 
this phenomenon could influence the equation of state at high densities achieved in heavy 
ion collisions or neutron stars. 

Below, we consider in detail the solutions of pion dispersion equation 

- k"^ -ml-U{iu,k,pF) = (1) 

in the complex plane of the pion frequency u. In (1) H is the pion polarization operator 
(pion self-energy) in the matter. The consideration in the complex cj-plane allows us to 
obtain additional information about well-known solutions. 
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The main content is as follows. We consider excitations with quantum numbers 0~ 
in symmetrical nuclear matter. At the equilibrium density there are three branches of 
solutions of (1). When the density is larger than critical one (p > pc) the fourth branch 
appears on the physical sheet. 

The equation (1) has logarithmic cuts on the physical sheet of the complex cu-plane, 
determined by the structure of the polarization operator H. Different branches of solutions 
are considered on physical and unphysical sheets of the complex cj-plane. With the change 
of k they move from physical to unphysical sheet (and backward) through the cuts. Below, 
to explain the obtained results, we deal with the case when isobar width in medium is 
equal to that in vacuum, i.e. 115 MeV. In this case all solutions of (1) and certain cuts 
move from the real (or imaginary) axis to the complex plane. This helps us to trace the 
^-dependence of the solutions. For certain important cases the influence of on the 
behaviour of uji{k) is studied. 

The well-known branches of solutions (1) are: 1) spin-isospin sound branch uJs{k)] 2) 
pion branch a;^(/c); 3) isobar branch UA{k). At < k < kf they are located on the physical 
sheet (the momentum kf is different for each branch) , then with further increase of k they 
move to unphysical sheet across the cuts. For these branches Iie{ujf) > everywhere on 
the physical sheet. 

Our investigations show that there is the whole set of solutions at unphysical sheets. 
With increasing density certain solutions come to the physical sheet. In such a way the 
fourth branch, udk), appears on the physical sheet at pp > 283 MeVQ. Below LJc{k) is 
referred as the condensate branch. At k = Ucik) is located at the same point as oj-^ik)] 
with increasing of k the branch goes onto the unphysical sheet. At some k = ki'it appears 
on the physical sheet and at k = k2 leaves it for the same unphysical sheet. The values of 
ki and ^2 depend on the density. At critical density, pc, corresponding to pp = 283 MeV 
(pc — l-2po) there is the equality ki = k2 = 1.8m.,r- It is Uc which obeys the inequality 
Recu^ < 0. All branches depend on the isobar width: at Fa = the branch ujc{k) is 
pure imaginary and uj'^ < 0. Recall that po is the equilibrium density of the matter, 
po = 2p^pq/3tt^, pfo = 268 MeV. 

The paper is organized as follows. In section 2 the particle-hole polarization operator 
n is considered in the complex ti;-plane. Then the branches of solutions of the pion 
dispersion equation (1) are presented on the physical and unphysical sheets. It is shown 
how the branches of solutions go across the logarithmic cuts in the complex plane. The 
condensate branch ujc is shown in details at different densities p and isobar width F^. 

2 Polarization operator and its singularities 

Here we write down the formulae for the polarization operator used. Our expressions, as 
is shown below, differ in some points from the well-known expressions [|l], 0] . Only S- and 
P-waves of the nNN (ttNA) interactions are taken into account. In this case 11 is the 

^For our values of the parameters, which are presented below. 
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sum of scalar, 115, sjid vector. Up, terms: 

U{uj,k) = Us{uj,k) + Up{uj,k) 



(2) 



The scalar polarization operator 115 is constructed using linear PCAC equation ob- 
tained by Gell-Mann-Oakes-Renner (GMOR) for the pion mass squared in the matter: 



,2 _ {NM\qq\NM){m^ + md) 



The value k = {NM\qq\NM) is the scalar quark condensate calculated in nuclear matter 
0; mu,md are masses of the current u- and d-quarks; /* is the pion decay constant in 
medium. Here we put / * = /tt = 92 MeV (more details can be found in 0). 

The scalar quark condensate n can be expanded in a power series with respect to p 

E,i 

K, = kq + p{N\qq\N) + terms with higher degrees of p, (4) 

where kq is the value of scalar quark condensate in vacuum, kq = —0.03 GeV'^; (A^|gg|A^) 
is the matrix element for scalar quark condensate in nucleon, (A^|gg|A^) ~ 8. In this place 
it is enough to keep the first two terms in (4), this correspond to the gas approximation 
for K M. Then we get from (3) 



m 



*2 2 {N\qq\N){mu + md) 

- ^n-P 2P • 



On the other hand, in the dispersion equation (1) m*^ is defined as 

mf = ml + n(cj, k = 0) . (6) 

Whilst Ilp{k = 0) = 0, we get from (5) and (6) the following form for 115: 

{N\qq\N){mu + md) 
Ils = -P 2^ . (7) 

The P-wave polarization operator Hp can be written following the papers |1|, 0, ^, [1^] 
as a sum of nucleon and isobar polarization operators 

Up = + Ha . (8) 

Here Hat (Ha) is equal to the sum of the nucleon-hole and isobar-hole loops without pion 
in the intermediate states H, H, H] : 

^ ^0 1 + (7a-7aa) n°,/fc^ l + (7A-7iV7v) 11%/^ 
iljv = Hat ^ ) = i^A ^ > (yj 



E . ^ ^ ^ 

Hat / 2 

E = 1- iNN-r^ - lAA-r^ + [iNNlAA " 7a 



k^ 
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The expressions for the nucleon-hole loop , 11%, is as follows: 



(27 



tNN 



6{p - Pf)6{pf - \ p + k\) _^ d{PF - p)6{\p + k\- pp) 



Here Ep = p^ /2m*. The expression for 11^ is analogous. 

The T[NB vertex T^^p^b with B labelling nucleon or A-isobar is 



ttNB 



tNB 



^0 

nNN 



9a 



9A 



(10) 

(11) 
(12) 



where x is nucleon 2-spinors, a is nucleon spin, turns the spin 3/2 into 1/2. In order 
to take into account the non-zero baryon size, the vertex r°^^ is multiplied by the form 
factor dsik) taken in the form (Ib = {I — m^/A|)/(l + A;^/A|). 

The constants 7ArAr,7A,7AA are 



7nn 



. 9A 



.7A 



Co9na ( /tt 

/A/Af 



9A 



,7aa 



Co9'aa ( /tt 



A/N 



9A 



where Co is the normalization factor (Co = ir'^ /{pfiti*)) and g's are the constants of the 
effective quasi-particle-quasi-hole interaction in nuclear matter |l|, ^. The calculation 
results are given for the following set of parameters |1|, |^, ^ : 

U/N ^2, = 0M7GeV, Aa = IGeV, qa = 1, f^ = 92MeV, (13) 

9'nn = 1-0, fl-^vA = 0-2, 9'aa = 0-8. 

When the parameters fA/N,9',A.B vary within the experimentally acceptable limits the 
dispersion equation solutions change quantitatively but not qualitatively. 

Later on, we can perform integration in (10) in two ways. 

1. To integrate separately the first and the second terms. In this way there appears 
the expression for 11^ as follows: 

/ « \ ^ 

Il%{u,k) = -4 



m 



^Ni(^,k) = — 



( 3^ 
1 f- 



) e[^N{u;,k) + ^N{-uj,k)]dUk) 
■uim* + kpF 



(14) 



k 47r2 



um* In 



^ {kpF^ - {ujm* - P/2y J / um* 



^iom* — kpF + k"^ 1'^ , 
- kpF - A;V2\\ 



+ 



2k'^ \ijjm* — kpp + k"^ /2 ^ 

at < < 2pF. At k > 2pF $Ar(co', k) is Migdal's function: 



(15) 



^Ni^,k) 



*3 



1 m 
4^ "P" 



In 



' a + 
a — b , 



ab 



(16) 
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where a = uj — (A;^/2m*), h = kpF/m*. Consider now which cuts has the polarization 
operator Il%{uj,k) (14)-(16) in the cu-plane. We can see that at < 2pF there are two 
cuts (define them as I and //). They are related to the first and second logarithms in 
(15). The cuts are situated within the intervals: 

T n ^ ^ TT ^ ^ ^Pp , (^^\ 
I : < UJ < , II : < < 1 • (17) 

m* 2m* m* 2m* m* 2m* 

Since 11^ is symmetrical under the replacement uj ^ —u, the cuts of ^N{—uJ,k) are 
placed symmetrically on the negative semiaxis. Thus 11^ has four cuts in the complex 
cu-plane, they are shown in Fig.l. 

2. The other way of integration in (10) gives a well-known expression of 11^ through 
Migdal's functions. To follow it, let us do a substitution in (10) O^p—pp) l — 9{pF—p) = 
1 — n{p). Then, after the integration, <I>Ar(a;, k) takes a well-known form (16) for all values 
of k [|l], ^ ^. Undoubtedly, the expression (14) is the same for both integration ways. 
Now has not four cuts in cj-plane but two (overlapping). It can be seen from the 
expressions (10), (15), that the cut / is the sum of two overlapping cuts. The search for 
dispersion equation solutions becomes more convenient and obvious when we work with 
one cut but not with two overlapping ones. It is difficult to follow the solution in the case 
2, therefore we use equations (14)-(16) for Il^^r. 

Now let us turn to the polarization operator 11^, which is the isobar-nucleon-hole loop. 
Since the isobar Fermi surface is absent at the nuclear densities and isobar momentum is 
unrestricted, there is no problem discussed above and n^(a;, k) reads: 

2 



^A{i^,k) + ^A{-UJ,k) 



di{k) . {U 



The function $A(t^, k) is expressed through Migdal's functions (16) with a = u)—{k'^ /2m*) — 
Am, h = kpp/m* . The mass difference. Am = m/\ — m, is the following: Re(Am) = 292 
MeV and Im(Am) = — rA/2. The cuts of ^) shown in Fig.l. At > the cut 

is in the interval 

^ + A^_^ < , < ^ + Am + ^. (19) 
2m* m* 2m* m* 

The cut is shifted into the complex plane in the value — zrA/2. 

In this paper Landau equation |ll|] is used for the nucleon effective mass 

= i^,o / 2w • (20) 
1 + {2mpF/TX^)fi 

Unknown parameter /i is fixed by the condition m*{pp = ppo) = 0.8m. 



3 Solutions of the dispersion equation 

In this section the solutions of the dispersion equation (1) are presented. The solution 
branch Uc emerges on the physical sheet of the complex w-plane at pi? = 283 MeV for the 
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parameter values (13). The figures for tlie zero sound brancli uJs{k), pion brancli uj.,r{k) 
and isobar brancli uj/^{k) are presented for Fa = 115 MeV at pp =268 and 290 MeV (i.e. 
at the equilibrium density and at density slightly larger than critical one). Some special 
cases, with the other values of pp and Fa, are considered as well. 

Branch u!s{k). (Fig.2a) The branch uJs{k) is shown for pp = 268 and 290 MeV 
(curves 1 and 2 correspondingly). The branch begins at uJs{k = 0) = 0, then while k 
increases, moves practically along the real axis. At kf = 0.430m7r for pi? = 290 MeV (at 
kf = 0.436m^ for pp = 268 MeV) goes under the cut // (17), this corresponds to the 
decay of into real nucleon and the hole. 

Branch ujA{k)- (Fig.2b) The isobar branch uji^{k) begins at = Am at = and 
ends on the isobar cut (19) at kj = 5.1m^ for pp = 268 MeV {kf = 4.8m^ for pp = 290 
MeV). 

Branch uj.^{k). (Fig. 2c) The pion branch starts at /c = in cJtt = m* (see (5), (6)). 
The beginning of uj^^lk = 0) is shifted to the smaller then values when GMOR P 
is used to determine m*. The pion branch ends on physical sheet under the isobar cut, 
this corresponds to the decay of pion into isobar and nucleon hole. It takes place at 
kf = 3.5m^ for pp = 268 MeV {kf = 3.8m^, pp = 290 MeV). 

Branch ujc{k). (Fig.2d, 3a,b) While the density increases one more branch of solu- 
tions, u!c{k), appears on the physical sheet. It emerges at pp > 283 MeV when parameter 
values (13) are used. In Fig. 2d the pion branch u!T,{k) and condensate branch u!c{k) are 
presented at pp = 290 MeV. The dashed piece of u!c{k) belongs to the upper unphysical 
sheet of the logarithmic cut / (17) (Fig.l). The branch u!c{k) starts at /c = at the same 
point as u!.^{k) and moves onto the unphysical sheet; at k = ki = 1.3m^ the branch goes 
down to the physical sheet and at = ^2 = 2.3m7r moves back to the same unphysical 
sheet. In the momentum interval {ki,k2) one can follow over all the branches shown in 
Fig. 2 and check that Uc does not belong to any branches studied before {us,Ut,,uj^). 

The branch Uc depends on the isobar width Fa. Decreasing the isobar width we see 
that the isobar cuts move to the real axis and ujg, uj^ and cja have smaller imaginary 
parts. When Fa = the branches cu^, uj.„ and wa are real. On the contrary, ujc{k) moves 
to the imaginary axis with decreasing Fa (Fig. 3a). At Fa = we have pure imaginary 
solutions on the physical sheet: cu^ < 0. 

In Fig. 3b the branch uJc{k) is shown at the different densities: pp = 280, 290, 300, 360 
MeV. When pp = 280 MeV the whole branch (curve 1) is located on unphysical sheet. 
At Pp = 283 MeV the branch touches the real axis (not shown). For pp > 283 Mev the 
part of u!c{k) in the interval {ki, /C2) is placed on the physical sheet. At the critical density 
p = rhoc and Fa = uJc{k) touches the real axis at the point ujc{k) = 0. 

It was shown in paper [0] that the appearance of Uc on the physical sheet results not 
only in pion condensation but in restoration of chiral symmetry in the nuclear matter at 
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critical density as well. 



4 Conclusion 

In the paper the solutions of pion dispersion equation are considered in details in the com- 
plex cu-plane. It is shown that, besides the well-known solutions with quantum numbers 
0~ (zero spin-isospin sound, pion and isobar waves), there exists the fourth branch u!c{k). 
It is the branch which obeys the condition o;^ < 0, therefore it is responsible for instability 
of the ground state. Such instability can indicate the beginning of 'pion condensation'. 
We demonstrate that at the density less than critical one, p < the branch U}c{k) is 
situated on unphysical sheet and at p > Pc it comes on physical one. 
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5 Figure captions 



Fig.l. The cuts on the physical sheet of the complex a;-plane of polarization operators 
n^, n^, corresponding to equations (14), (17), (18), (19). The cuts are presented at 

Pf = 290 MeV, k = m^. 

Fig. 2. Branches of solutions of (1) in the complex cj-plane. Curves 1 and 2 stand for pp = 
268, 290 MeV, correspondingly. The dashed pieces of curves are situated on unphysical 
sheets, a) Zero spin-isospin sound branch uigik). The curves are presented up to 1.6m^. b) 
The isobar branch loa- The horizontal dashed hue is a logarithmic cut (19) for = 290 
MeV at momentum k when uiA{k) is on the cut. c) The pion branch cj^r- The horizontal 
dashed line is a logarithmic cut (19) for pp = 290 MeV at momentum k when cj^ is on 
the cut. d) The total picture at pp = 290 MeV for pion branch a;^(fc) and condensate 
branch ci;c(k). 

Fig. 3. Condensate branch cUc in the complex a;-plane. Here the dashed pieces of branches 
are on the physical sheet, but solid lines belong to unphysical sheet, a) The branch ujc is 
presented at pp = 290 MeV for different values of isobar width: Fa = 0, 10, 50, 115 MeV 
(curves 1,2,3,4, correspondingly); cudk = 0) = 0.744m^. b) The branch Uc is presented 
at Fa = 115 MeV for different values of Fermi momenta pp = 280, 290, 300, 360 MeV 
(curves 1,2,3,4, correspondingly). For pp— 300 and 360 MeV only the part of the branch, 
which is placed on the physical sheet, is shown (the curves 3 and 4). The whole branch 
for Pp — 280 MeV (curve 1) is on unphysical sheet. 
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